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We studied magneto-optical resonances caused by excited-state level crossings in a nonzero mag- 
netic field. Experimental measurements were performed on the transitions of the D2 line of rubidium. 
These measured signals were described by a theoretical model that takes into account all neighboring 
hyperfine transitions, the mixing of magnetic sublevels in an external magnetic field, the coherence 
properties of the exciting laser radiation, and the Doppler effect. Good agreement between the 
experimental measurements and the theoretical model could be achieved over a wide range of laser 
power densities. We further showed that the contrasts of the level-crossing peaks can be sensitive to 
changes in the frequency of the exciting laser radiation as small as several tens of megahertz when 
the hyperfine splitting of the exciting state is larger than the Doppler broadening. 

PACS numbers: 32.60.+i,32.80.Xx 



I. INTRODUCTION 

Level-crossing spectroscopy has long been used to 
study lifetimes of atomic states (using zero-field reso- 
nances) or atomic constants, such as the magnetic mo- 
ments and fine and hyperfine constants (using nonzero- 
field resonances) [U [2] . The technique most often takes 
advantage of resonances in plots of the laser-induced flu- 
orescence (L1F) in a particular direction with a given 
polarization as a function of the magnetic field. The res- 
onances are related to the type of coherent excitation 
of magnetic sublevels that becomes possible when some 
of them, whose z-components mj? of the total angular 
momentum F differ by Amp — q, become degenerate at 
particular magnetic field values [3] . For linearly polarized 
excitation q = ±2. Such a degeneracy always occurs at 
zero magnetic field where all magnetic sublevels belong- 
ing to a particular hyperfine level F have the same energy. 
This case of zero-field level crossing is known as the Hanlc 
effect, first observed by Hanle himself [J. However, as 
can be seen in Fig. [TJ it also happens at certain nonzero 
magnetic field values that some magnetic sublevels from 
different hyperfine F states can cross. If the requirement 
for coherently excited magnetic sublevels with a certain 
Am is fulfilled, one speaks of nonzero-magnetic-field level 
crossings. The coherent evolution of such systems can be 
described with the optical Bloch equations (OBEs) for 
the density matrix. However, in order to describe accu- 
rately real systems, it is necessary to take into account 
all neighboring hyperfine transitions, the magnetic-field- 
induced mixing of magnetic sublevels of identical m that 
belong to different hyperfine levels, the Doppler profile, 
and the coherence properties of the radiation. Models 
with these characteristics have been developed over the 
years to describe zero-field resonances in the ground and 
excited states with great precision [5]. In this work, we 
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show that nonzero level-crossing signals in magnetic fields 
can be described by a theoretical model over a wide range 
of magnetic fields to nearly experimental accuracy. More- 
over, the model succeeds also at laser power densities for 
which the excitation is nonlinear and where the effects 
of optical pumping can be noted. We also show that, 
by carefully selecting and controlling the laser frequency 
with a precision of several tens of megahertz, it is possible 
to obtain resonances with greater contrast, which can be 
useful in applications such as measurements of magnetic 
fields in the range of tens of Gauss or determination of 
hyperfine constants in states for which they are not yet 
known. The ability to describe level-crossing signals pre- 
cisely can be useful for determining atomic constants es- 
pecially in situations where the large number of crossing 
points washes out individual resonances. 

The first theory of level-crossing signals was given by 
Breit in 1933 [S], and the first application was to mea- 
sure the fine structure splitting between the helium P 
states [7J. These measurements were described in terms 
of Breit's formalism by Franken in 1961 [5]. Since then, 
these signals were used extensively for a time to make 
measurements of the fine and hyperfine constants in 
atoms. For example, the technique was used to obtain 
hyperfine constants in rubidium [5] and cesium [10] (see 
Ref. p] for a review of many results). Theoretical mod- 
els of the ground-state Hanle effect were used by Picque 
in 1978 [TT]. Over time, these models became more and 
more sophisticated as different effects were included [H] . 
Precise analytical models are also possible [13], but only 
for lower laser power in the linear regime. 

The present work revisits an earlier study published in 
2003 [14 , which offered only a model that was limited 
to the cycling transitions in the limit of weak excitation 
and thus could provide only a qualitative description of 
the experimental signals. We now show that the signals 
can be described very precisely even in the case of strong, 
nonlinear excitation with a model based on the optical 
Bloch equations and valid also for nonlinear excitation. 
Our model takes into account possible contributions to 
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the transition probabilities from all neighboring hyperfine 
transitions, the effects of Doppler broadening, the split- 
ting of the hyperfine levels in the magnetic field, and the 
coherence properties of the exciting laser radiation |15j . 
This model had been widely applied to zero-field reso- 
nances in the ground state [5] and achieved good agree- 
ment between experimentally measured and calculated 
curves. The experimental parameters, in particular the 
laser frequency, were carefully controlled during the mea- 
surements in order to allow precise comparison with the- 
ory. In addition, experimental measurements and theo- 
retical calculations were used to investigate the influence 
of the laser frequency on the relative contrasts of the 
level-crossing peaks at nonzero field values. 



II. EXPERIMENT 

Figure [T] shows the relative energies of the excited state 
magnetic sublevels as a function of the magnetic field. 
Each curve corresponds to a particular value of the 
projection of the total angular momentum F on the z- 
axis. As we used linearly polarized exciting radiation, co- 
herences can be formed around the crossing points with 
Am = ±2, which are circled and labeled by small Greek 
letters. During the experiments, the laser frequency de- 
tuning is determined as the frequency difference between 
a particular ground-state hyperfine level and the excited- 
state fine structure level. 

In this experiment natural atomic rubidium, confined 
in a Pyrex cell with optical quality windows (25 mm long 
and 25 mm in diameter), was placed in the middle of a 
three-axis set of Hclmholtz coils. Two pairs of coils were 
used to compensate the Earth's magnetic field, while the 
magnetic field was scanned along the third (z) axis. As 
the magnetic field was scanned through a triangular pat- 
tern with a frequency of 0.02 Hz, fluorescence spectra 
were acquired. The laser wavelength was determined by 
means of a saturated absorption spectroscopy setup in 
conjunction with a WS-7 wavemeter from HighFinesse. 
It was monitored during the scan with the wavemeter, 
and adjustments were made if necessary. Using a bipolar 
Kepco BOP-50-8-M or an Agilent N5770A power sup- 
ply, magnetic fields of up to 120 G could be achieved. 
Laser radiation from an external cavity diode laser passed 
through a chopper and entered the cell with its propa- 
gation vector and electric field vector both perpendic- 
ular to the scanning magnetic field (see Fig. [2]). The 
temperature of the laser box and the diode were stabi- 
lized by Thorlabs TED200 temperature controllers and 
the current was controlled by a Thorlabs LDC205B cur- 
rent controller. The diameter of the beam was 1.6 mm as 
measured by a Thorlabs BP104-VIS beam profiler. The 
beam width was defined as the full width at half max- 
imum of the Gaussian intensity profile. By means of a 
polarization rotator followed by a linear polarizer, laser 
power values from 20 pjW to 320 pW could be achieved, 
which translated into laser power densities of 1 mW/cm 2 



to 16 mW/cm 2 . The LIF of two mutually perpendicu- 
lar components (one parallel and the other perpendicu- 
lar to the exciting electric field vector) was detected by 
two Thorlabs FDS100 photodiodes located behind a po- 
larizing beam-splitter. The signals were amplified by a 
transimpedance amplifier followed by a lock-in amplifier 
and recorded using an Agilent DSO5014A oscilloscope, 
which also averaged the signals. 



III. THEORETICAL MODEL 

We describe the atomic system via its quantum den- 
sity matrix p, which is written in the basis of Zeeman 
sublevels for the hyperfine structure of the Di transition 
of atomic Rubidium: Fi, rriFi), where Fj denotes the 
quantum number of the total atomic angular momentum 
in either the ground (i — g) or the excited (i = e) state, 
rriFi refers to the respective magnetic quantum number, 
and £i represents all other quantum numbers that remain 
unchanged within the Di line. The time evolution of the 
density matrix p is governed by the optical Bloch equa- 
tions [TBI: 



lh di = 



H,p 



ihRp, 



(1) 



where H denotes the full Hamilton operator of the sys- 
tem and R is the relaxation operator. The full Hamilto- 
nian can be expressed in terms of the unperturbed atomic 
Hamiltonian H determined by the internal dynamics of 
the atom, a term Hb that describes the interaction with 
the external magnetic field, and a dipole interaction term 
y = -d-E(t): 



H = H + H E 



V. 



(2) 



As indicated, the interaction with the electromagnetic 
field is treated in the dipole approximation [17] . The 
magnetic interaction Hamiltonian can be written as 



Pb , 



jjj + <7/I) ■ B, 



(3) 



where ps is the Bohr magneton, J and I are the total 
electronic angular momentum and spin of the atomic nu- 
cleus and <7j, gj are the respective Lande factors. The 
interaction Hamiltonian ^ consists of the interaction 
matrices for fixed projection mp = mj + mj of the an- 
gular momenta on the quantization axis, which can be 
written in terms of Wigner 3 j symbols |17j . Solving the 
eigenvalue problem for these matrices yields the energy 
structure shown in Fig. [TJ where J and I make up the 
total atomic angular momentum F = J + I. The com- 
ponents (cfe) of the eigenvectors are used as mixing co- 
efficients for the chosen basis. The mixed atomic states 
have to be rewritten in the presence of the magnetic field 
as 



\Z,F,m) 



^2c k \£,F k ,m). 



(4) 
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FIG. 2. (Color online) Schematic drawing of the experimen- 
tal setup. 



The density matrix p can be divided into four parts: 
p gg and p ee are quadrants aligned along the main di- 
agonal and are called Zeeman coherences, p ge and p eg 
are called optical coherences. In order to describe the 
atomic fluorescence one has to know the elements p ee - 
The atomic dipole interaction with a classically oscillat- 
ing electric field characterized by a stochastic phase that 
gives raise to some spectral width Aw of the frequency 
is considered. The rotating wave approximation [18] and 
the decorrelation of and averaging over the stochastic 
phase [IH] allows adiabatic elimination of the optical co- 
herences |15j , which results in rate equations for the Zee- 



man coherences: 
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^Pe.eS)- 
i^e z e J Pe i e J ~ (T + l)Pe z e 3 ■ 



(5b) 

The following terms are used in Eq. ([5]): describes 
the interaction strength between the atom and the laser 
radiation and is defined below, dij is the dipole transi- 
tion matrix element that can be obtained from the re- 
duced matrix element by means of the Wigner-Eckart 
theorem |17j . Uij is the energy difference between lev- 
els \i) and \j), Tg*p. describes coherence transfer to the 
ground state via spontaneous emission, 7 is the transit 
relaxation rate at which the atoms leave the interaction 
region, T is the rate of the spontaneous transitions, and 
A describes rate at which "fresh" atoms enter the inter- 
action region. We assume that atoms entering the in- 
teraction region are completely depolarized and that the 
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atomic equilibrium density outside the interaction region 
is normalized to unity; thus, A = 7. 

The interaction strength Sjj is given by 

*.. = ]£>»!! 1 (6) 

where £q is the amplitude of the oscillating electric field, 
Q denotes the central frequency of the laser radiation, Aw 
is the laser linewidth, and k s v is the Doppler shift of the 
atomic transition for an atom moving with velocity v. In 
the numerical simulations we define the Rabi frequency 
in the following way: 

\e a \-\\dj\\ m 

= h ' 

where \ \dj\ \ is the reduced dipole element of the D2 tran- 
sition whose value can be found in [17J . The Rabi fre- 
quency is used as a parameter that corresponds to the 
experimentally measurable excitation laser power den- 
sity: 

/ = k Rabl n 2 R . (8) 

The proportionality coefficient kji a {,i can be estimated 
from Eq. Q. However, the precise value must be deter- 
mined from the best fit to the experimental data because 
the laser power density is not constant over the beam 
profile. 

In this study we numerically solved equations §5j for 

steady state excitation conditions ( dP g t S3 — —^f 1 - = o) 

to obtain the Zeeman coherences p gg and p ee . From this 
point it was straightforward to obtain the fluorescence 
for some particular polarization component defined by e: 

I fl (e)=I £ 4?>dW Pejeh , (9) 

where d[° b ^ are the dipole transition matrix elements for 
the chosen observation component. 

In order to take into account the classical movement 
of atoms with velocity v we performed a numerical in- 
tegration of the fluorescence signal over the frequency 
distribution in the Doppler profile. 

IV. RESULTS AND DISCUSSION 

As a first test of how well the model can describe non- 
zero level-crossing signals, we used the model to repro- 
duce the experimental results of non-zero level crossing 
signals for the F g = 2 — > F e = 1, 2, 3 transitions of 85 Rb. 
The results are shown in Fig. [3j Measurements were car- 
ried out for a range of laser power densities I from 1 
mW/cm 2 to 16 mW/cm 2 . The Rabi frequency £Ir was 
assumed to vary as in Eq. pi, where the same value for 
kRabi was used for all calculations in this paper. This 



value was chosen to provide the best overall agreement 
while being consistent with a rough estimation based on 
theoretical considerations (see Eq. Q). The other pa- 
rameters that were optimized in a similar way were the 
laser linewidth Aw and the proportionality constant be- 
tween the laser beam width and the transit relaxation 
time. In each case only one constant value was used in 
all the theoretical calculations presented in this work. 
The agreement between theory and experiment was ex- 
cellent at low laser power densities and still quite good 
at / = 16 mW/cm 2 . At higher laser power densities the 
model becomes less accurate, because at higher power it 
is no longer sufficient to describe the transit relaxation 
by a single rate constant [20] as this model does. Never- 
theless, the model seems adequate to predicting signals 
for level-crossing spectroscopy with Am = ±2. 

Next we turned our attention to the influence of the 
laser detuning on the level-crossing spectra. Figure [4] 
presents the results of theoretical calculations that show 
how the relative contrasts of the nonzero-field level- 
crossing peaks can be influenced by the tuning of the 
exciting laser radiation. The detunings of the different 
curves plotted are measured with respect to the energy 
difference between the indicated ground-state hypcrfine 
level and the excited-state fine structure level. First of 
all, it was apparent that significant variations in the rela- 
tive contrast of the level-crossing peaks could be achieved 
only in 87 Rb [Fig. |4^a,b)]. The contrast of the peak at 
level crossing a could be increased by up to a factor of 
four when exciting from the F g =2 ground state. When 
exciting from the F g =l ground state, the contrast of the 
peak that corresponds to level crossing j3 can be increased 
by almost a factor of two. In the case of 85 Rb [Fig.gc,d)] 
Doppler broadening washes out any possible effects, be- 
cause the hyperfme splitting of the excited state in this 
isotope is small relative to the Doppler width. The en- 
ergy difference between the F = 1 and F = 4 states of 
85 Rb is slightly more than 200 MHz, whereas the full 
width at half maximum of the Doppler profile is around 
500 MHz at room temperature. Although the change in 
laser frequency was smaller in absolute terms for the cal- 
culations with 85 Rb than for the calculations with 87 Rb, 
the change in detuning as a fraction of hyperfme split- 
ting was roughly equal. Detuning the laser by a larger 
amount for 85 Rb would result in significantly lower sig- 
nals without affecting peak contrast significantly. 

Next we compare experimental level-crossing curves 
obtained at different values of the laser detuning with 
the calculated curves for those same detuning values to 
see if the predicted changes in contrast can be observed. 
Figure [5] shows results for 85 Rb. Figure [HJa) shows the 
results of an experiment in which the laser was detuned 
from the the exact energy difference between the ground 
state with F g = 2 and the 5 2 P 3/2 state by 65 MHz. The 
results shown in Fig. [5jb) were obtained with a laser 
detuning of —44 MHz. Similarly, Fig. |6](a) corresponds 
to a measurement with the laser detuning between the 
ground state of 85 Rb with F g — 3 and the 5 2 P 3 / 2 state 
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FIG. 3. (Color online) Intensity difference (I± — versus magnetic field for the F g = 2 — > F e = 1, 2, 3 transitions of 85 Rb. 
The different panels correspond to different laser power densities. Markers represent the results of experimental measurements, 
while the curves represent the results of theoretical calculations. 



of —63 MHz, while Fig. [6jb) corresponds to an experi- 
ment with a laser detuning of 22 MHz. As expected, the 
relative contrasts of the level-crossing peaks are not very 
sensitive to detuning. We note that the theoretical model 
describes the experimental curves quite well. These tran- 
sitions also illustrate the necessity of including magnetic 
sublevel mixing in the theoretical model. For example, 
the resonance at position e (see Fig. [T]) involves a cross- 
ing of the sublevels labeled by F e = 3, rap = —1 and 
F e — 4,mp = —3. Notwithstanding the selection rules 
for B = (AF — 0, ±1), this level crossing can produce a 
resonance even when the excitation takes place from the 
ground state with F g = 2. The reason is that F ceases to 
be a good quantum number when the magnetic field is 
nonzero, and, at high values of the magnetic field, there 
is a sufficiently high probability of exciting the state la- 
beled by F e = 4,mj? = —3 from F g = 2. However, the 
state with maximum projection of angular momentum 
[F e = 4,mp = 4) is not mixed. Thus, it is not possible 
to excite it from F g = 2, but only from F g = 3. These 
considerations are borne out by the theoretical and ex- 
perimental results shown in Fig. [5]and Fig. [6] A theoret- 
ical model, such as the one in [j"4] [ that did not take into 
account magnetic sublevel mixing at high magnetic fields 



would not have been able to reproduce the resonance at 
position e for an excitation from F g = 2 [14] . 

Figure[7]and Fig.|8]show results for 87 Rb for excitation 
from the F g — 1 and F g = 2 ground state levels, respec- 
tively. In this case, the detuning can dramatically affect 
the shape of the signals. As can be seen by comparing 
Fig.[7]Ja) and Fig.[7][b), the contrast of the peak that cor- 
responds to level crossing /3 increased by more than a fac- 
tor of two when the laser was detuned by 36 MHz as com- 
pared to when it was detuned by —190 MHz. Similarly, 
when the excitation took place from the ground state 
with F g — 2, the peak that corresponds level crossing a 
was not visible when the laser was detuned by 36 MHz 
[Fig. [8]ja)], but appeared when the laser was detuned by 
-190 MHz [Fig. [§]»]. The contrast of the peak at level 
crossing 7 also increases by a factor of two depending on 
if the laser is tuned to a value that closely corresponds 
to the energy difference between the ground state and 
its level crossing [Fig. |8jb)] or far away [Fig. [8^a)]. Note 
that "close" here does not mean in magnetic field val- 
ues, but in energy difference (see Fig. [I]) Again, these 
peaks illustrate the power of using a detailed model, as 
the resonance labeled /? would not have appeared in the 
theoretical calculation when exciting from F g = 1 with- 
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FIG. 4. (Color online) Theoretical calculations of I± — Jy versus magnetic field for (a) the F g = 2 — > F e = 1,2,3 transition 
of 87 Rb, (b) the F g = 1 — > F e = 0,1,2 transition of 87 Rb, (c) the F g = 2 — ► F e = 1,2,3 transition of 85 Rb, and (d) the 
F s = 3 — > F e = 2, 3, 4 transition of 85 Rb. The different curves in each figure correspond to different laser frequencies. The 
peaks are labeled by small Greek letters to indicate which crossing points gives rise to them (see Fig. [TJ . 




FIG. 5. (Color online) Intensity difference (Jj_ — versus magnetic field for the F g — 2 — > F e — 1,2,3 transitions of 85 Rb 
with the laser detuned from the 5 2 S'i/2(f g = 2) — > 5 2 P3/2 transition by (a) 65 MHz and (b) —44 MHz. Markers show the 
results of an experiment with the laser power density 1 = 2 mW/cm 2 , while the solid line shows the results of a calculation 
with Rabi frequency Q,r = 4 MHz. 
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FIG. 6. (Color online) Intensity difference (I± — versus magnetic field for the F g = 3 — > F e = 2,3,4 transitions of 85 Rb 
with the laser detuned from the 5 2 Si/2(F g — 3) — > 5 2 P3/2 transition by (a) —63 MHz and by (b) 22 MHz. Markers show the 
results of an experiment with the laser power density 1 = 2 mW/cm 2 , while the solid line shows the results of a calculation 
with Rabi frequency Q.r = 4 MHz. 



out taking into account sublevel mixing. 

Figure [9] shows the intensity difference as a function of 
laser detuning when the magnetic field is fixed to a value 
of B = 57 G and the excitation takes place from the 
ground state hyperfine level with F g = 1. One can clearly 
see that there is an optimal detuning for maximizing the 
intensity difference between the two polarization compo- 
nents. The curve has a half width at half maximum of 
around 250 MHz, comparable to the Doppler width of 
rubidium at room temperature. The black squares con- 
nected by a line give the results of a theoretical calcula- 
tion for a Rabi frequency of fin = 4 MHz. The filled cir- 
cles represent experimentally measured values for a laser 
power density of 2 mW/cm 2 . The agreement is quite 
good. In order to optimize the experimental conditions 
under which the maximum contrast for a level-crossing 
peak could be observed, it can be useful to generate from 
the results of the calculations a three-dimensional plot of 
the intensity difference of the two orthogonally polarized 
components of the LIF as a function of laser detuning 
and magnetic field. Such a plot is shown in Fig. [10] and 
from it one can easily determine the optimal laser de- 
tuning for maximizing the contrast of each level-crossing 
peak. 

As is well-known, the level-crossing resonances are 
manifestations of coherences in the excited state manifold 
of the atoms under study. These coherences can be visu- 
alized by plotting surfaces that represent the probability 
of finding the angular momentum of an atom pointing 
in a particular direction in space [STJ [22] . Such surfaces 
can be generated from the density matrices by summing 
over all F states (see [53]), which yields the spatial distri- 
bution of J. Such surfaces have been included in Fig. [7] 
as an illustration of the information that can be gleaned 
from the theoretical model. The axes I x and I y refer 
to the states that give rise to the respectively polarized 
fluorescence when they decay. When two or more sub- 
levels in the excited state that could be excited coherently 
(Am = ±2) are degenerate, the excited state becomes 
aligned, which manifests itself as a nonuniform angular 
momentum distribution in the xy-plane. In other words, 
an aligned state gives rise to an angular momentum dis- 
tribution that has the z-axis as a second-order symmetry 
axis. In the absence of coherence the angular momentum 
spatial distribution is axially symmetric with respect to 
the z-axis. Aligned states occur when all sublevels cross 



at zero magnetic field and at the level-crossing points 
that can be excited coherently, such as a and /? in Fig. [7] 
Far from the level-crossing points the angular momentum 
distributions become symmetric. The degree of asymme- 
try in the angular momentum distribution also reflects 
the strength of the coherence and the contrast of the re- 
spective resonance. 

V. CONCLUSION 

The results presented in this work demonstrate two 
things: (1) the shapes of Am = ±2 level-crossing sig- 
nals can be theoretically modeled with good accuracy 
over a broad range of magnetic field values and over at 
least an order of magnitude in laser power density; and 
(2) the detuning of the laser can dramatically influence 
the shape of the level-crossing curve, in particular the 
contrast of the peaks, when the hyperfine splitting of 
the excited state exceeds the Doppler broadening. This 
sensitivity to the laser detuning also provided a partic- 
ularly stringent test of the theoretical model, as well as 
being interesting in their own right for the optimization 
of level-crossing studies. Furthermore, the density matri- 
ces computed with the theoretical model can shed light 
on the coherent processes associated with the excitation 
of the atoms, which is useful for computing fluorescence 
intensities for arbitrary directions and polarizations. Al- 
though the atomic constants involved in the transitions of 
the Z?2 hue of rubidium are well known, a precise model 
of level-crossing signals can be necessary for extracting 
atomic constants in situations where the large number 
of level crossings washes out individual peaks or where 
the hyperfine splitting is small [23] . The utility of such a 
model for the hyperfine constants A for the 7,9, IO-D5/2 
states of cesium from electric field level crossings was 
demonstrated, for example, in |25j . 
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FIG. 7. (Color online) Intensity difference (I± — Iu) versus magnetic field for the F g = 1 — > F e = 0, 1, 2 transitions of 87 Rb 
with the laser detuned from the 5 2 Si/2(F g = 1) — > & 2 Pz/2 transition by (a) 36 MHz and by (b) —190 MHz. Markers show the 
results of an experiment with the laser power density 1 = 4 mW/cm 2 , while the solid line shows the results of a calculation 
with Rabi frequency Qr = 6 MHz. 
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FIG. 8. (Color online) Intensity difference (I± — Iu) versus magnetic field for the F g = 2 — > F e = 1, 2, 3 transitions of 87 Rb 
with the laser detuned from the 5 2 Si/2(F g = 2) — > 5 2 Ps/2 transition by (a) 36 MHz and (b) by —190 MHz. Markers show the 
results of an experiment with the laser power density 1 = 2 mW/cm 2 , while the solid line shows the results of a calculation 
with Rabi frequency J1_r = 4 MHz. 
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laser detuning for fixed magnetic field (B = 57 G) value for 
the transition from the ground state of 87 Rb with F g = 1. 
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